Using Gravitational Lensing to study HI clouds at high redshift 
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ABSTRACT 

We investigate the possibility of detecting HI emission from gravitationally lensed 
HI clouds (akin to damped Lyman-a clouds) at high redshift by carrying out deep radio 
observations in the fields of known cluster lenses. Such observations will be possible 
O ■ with present radio telescopes only if the lens substantially magnifies the flux of the HI 

emission. While at present this holds the only possibility of detecting the HI emission 
from such clouds, it has the disadvantage of being restricted to clouds that lie very 
close to the caustics of the lens. We find that observations at a detection threshold of 
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50 uJy at 320 MHz (possible with the GMRT) have a greater than 20% probability of 
detecting an HI cloud in the field of a cluster, provided the clouds have HI masses in 
the range 5 x 10 8 M < Mm < 2.5 x 10 10 M Q . The probability of detecting a cloud in- 
creases if they have larger HI masses, except in the cases where the number of HI clouds 
in the cluster field becomes very small. The probability of a detection at 610 MHz and 
233 MHz is comparable to that at 320 MHz, though a definitive statement is difficult 
owing to uncertainties in the HI content at the redshifts corresponding to these fre- 
quencies. Observations at a detection threshold of 2 ^Jy (possible in the future with 
the SKA) are expected to detect a few HI clouds in the field of every cluster provided 
the clouds have HI masses in the range 2 x 10 7 M < Mm < 10 9 M . Even if such 
observations do not result in the detection of HI clouds, they will be able to put useful 
constraints on the HI content of the clouds. 



subject headings: cosmology: gravitational lensing; theory: galaxies; large-scale structure of 
the universe; radio lines: HI 
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1. Introduction 

Quasar spectra reveal a variety of absorption features superposed on the continuum emission, 
indicating the existence of HI clouds with a wide range of neutral hydrogen column densities 
distributed at different redshift. Among these the damped Lyman-a clouds (DLA) which have HI 
column densities in the range 2 x 10 20 < Ahi(cui~ 2 ) < 10 22 are the main repository of neutral 
hydrogen at high redshift (z ~ 3). It has been speculated that these systems are the progenitor of 
present day spiral galaxies (for details see Wolfe 1995 and references therein). 

In addition to the column densities, absorption studies have determined the detailed velocity 
structure of the DLAs (Prochaska & Wolfe 1998), and different models have been proposed for 
DLAs in order to explain these observations. The DLAs have been modeled as thick rotating disks 
with rotation velocities in the range 200-300 kms -1 (Prochaska & Wolfe 1998). Another model 
proposes that DLAs are made up of gaseous protogalactic clumps undergoing in-fall (Haehnelt, 
Steinmetz k Rauch 1997). 

Despite very detailed absorption studies, the exact nature of the DLAs is still not fully under- 
stood. A part of the uncertainty arises from the fact that absorption studies give information along 
a single line of sight through a DLA, and this cannot be used to determine the total HI content or 
the physical extent of these systems. In a few cases there exist multiple lensed images of a quasar 
from which we can infer that the angular size of the DLAs is ~ 3" (Smette 1995 and reference 
therein) . Observations of DLAs in HI absorption and deep- imaging of the DLA fields in the optical 
band have also been used to infer the size of DLAs (Lane et al. 2000; Fynbo et al. 2000; Fynbo 
et al. 1999; Moller et al. 1998; Warren & Moller 1996; Moller & Warren 1993; Briggs et al. 1989; 
Briggs 1988), but such observations only give lower limits on the physical size of DLAs. 

In a few cases it has been possible to use deep optical imaging and spectroscopy to identify 
counterparts (in emission) of the systems that produce the damped Lyman-a lines. While most 
of these observations have been performed on low redshift DLA fields (see e.g. Brun et al. 1997 
and references therein) there also exist several cases of positive detection of DLAs at high redshift 
(Kulkarni et al. 2000; Djorgovski et al. 1996; Fynbo et al. 2000; Fynbo et al. 1999; Moller et al. 
1998; Warren & Moller 1996; Moller & Warren 1993). 

Here we consider the possibility of detecting HI emission from hydrogen clouds whose column 
densities are sufficiently high to qualify as DLAs, if they were located along the line of sight to 
quasars. This allows the properties of the HI clouds to be inferred from the observed properties 
of DLAs. These observations will also enable the total HI content of the DLAs to be determined. 
Unfortunately, the very small angular extent of the DLAs leads us to expect HI flux from these 
clouds to be below the detection threshold of existing radio telescopes. 

Optical observations of background galaxies (typically ^ 1") lensed by cluster of galaxies 
reveal arc-like images which could be as large as 20" (see e.g. Schneider et al. 1992). Detailed 
models have been constructed for these cluster lenses; and it is possible to predict the magnification 
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of the image of HI clouds if they were located at different positions behind the cluster. The lensing 
of HI clouds by clusters of galaxies holds the possibility of magnifying them sufficiently rendering 
the HI emission detectable by presently available radio telescopes. In this paper we investigate the 
possibility of detecting the redshifted 1420 MHz HI emission from lensed HI clouds by carrying out 
deep radio observations centered on known cluster lenses. 

Our investigation is largely motivated by the fact that the Giant Meterwave Radio Telescope 
(GMRT; G. Swarup et al. 1991) which has several bands in the frequency range 150-1420 MHz has 
recently started functioning. In this paper we consider frequency bands of width 16 MHz centered 
at 233, 320 and 620 MHz. These correspond to the redshifted 1420 MHz emission from HI at 
z ~ 1.3, 3.4 and 5.1 respectively. Most of the discussion in this paper refers to 320 MHz which we 
have used as the fiducial frequency. In addition we use a spatially flat Vt m = 1 model, and a Hubble 
parameter Hq = 100 hk.rn.s~ 1 . We shall use h = 0.5 whenever a numerical result is reported. 

The absorption studies of DLAs have yielded a maximum HI column density of 10 22 
(Lanzetta, Wolfe & Turnshek 1995). It is possible that clouds with higher column densities contain 
larger amount of dust that obscures the background quasar, in which case these clouds would not 
be detected in absorption (Fall & Pei 1995). This possibility has been taken into account while 
analyzing the possibility of detecting the HI emission by setting the maximum HI column density of 
the clouds to a value which is a few times larger than the maximum HI column densities observed 
in absorption. 

In § 2 we review the observed properties of HI clouds which are relevant for our analysis and 
describe two simple models for the HI cloud mass and flux distribution. In § 3 we briefly discuss a 
few relevant features of gravitational lensing by a cluster of galaxies and outline the method of our 
calculation. We present our results in §4 and in §5 we give our conclusions. 



We begin this section with a brief discussion of some of the observational features of DLAs 
which we use to determine the properties of the distribution of HI clouds. 

Absorption studies directly give the velocity dispersion AV of the HI along lines of sight 
through DLAs and these observations indicate rotational velocities in the range 200-250 kms -1 
(Prochaska & Wolfe 1998). This along with the HI column density can be used to calculate the 
specific intensity of the redshifted 1420 MHz emission from the HI along the line of sight through 
a HI cloud. This gives (for derivation see Appendix A): 



Integrating the specific intensity over the angular extent of the HI cloud gives the total flux from 
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Mass and flux distribution of HI clouds 




(1) 
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the cloud: 

F v <xl N m (J3) d 2 (3 , (2) 

where (3 (in arc-seconds) refers to different angular positions on the sky, and the integral is over 
the angular extent of the HI cloud. 

In a situation where there is gravitational lensing a source which would otherwise appear at 
(3 is seen at a different angular position 6 (see e.g. Schneider et al. 1992), and the flux from the 
lensed image of the HI cloud can be calculated by performing a similar integration over the angular 
extent of the image as in Eq. (2). 

Estimating the flux expected from a HI cloud (lensed or unlensed) requires knowledge about 
its angular extent and the variation in the column density across the HI cloud. As these facts are 
not available, we have adopted two simple models for our calculations in this paper. We describe 
these models below. 



2.1. Uniform disk (UD) 

In this model the HI clouds are assumed to be face on circular disks, all of the same physical 
radius r^LA and with uniform column density across the disk. We consider different values of t^la 
in the range 10 kpc to 30 kpc. The solid angle subtended by a HI cloud is AOdla = ^dla/^s ( z ) 
where d s (z) is the angular diameter distance to the HI cloud (see e.g. Peebles 1993). In this model 
the flux from an unlensed HI cloud is F v = /^ASIdla and we have, for v = 320 MHz, 

The distribution of the column densities of the HI clouds is assumed to be a power law F(Nm) = 
BN R * with values in the range 2 x 10 20 cm~ 2 < A H i < 10 22 cm~ 2 , where F(Nui) dNm gives the 
number of HI clouds per unit volume with column density in the range dNm- Observations indicate 
that the column density distribution changes with redshift implying that both <x and B change with 
z. The DLA population shows a marked increase in the HI column densities at higher redshifts 
(Lanzetta et al. 1995) and 75% of the DLAs in the redshift range 3 < z < 3.5 have column densities 
Nm > 10 21 cm -2 . This is consistent with a = 0.5. In our work we consider values of a in the range 
0.5 to 1.0. 

In this model the HI mass of a cloud is given by 

This relation permits us to interchangeably use either a column density distribution function or a 
HI mass distribution function /(Mhi) = where /(Mhi)g?Mhi is the number of HI clouds per 
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unit volume with HI mass in the range g?Mhi- The mass range corresponding to the column density- 
range 2 x 10 20 cm" 2 < A^hi < 10 22 cm -2 depends on the value of tdla> and for tola = 10 kpc 
we have Mnifmin] = 5 x 10 8 M Q and Mhi [max] = 2.5 x 10 10 M . At any redshift z the HI mass 
distribution function can be used to calculate the HI density 

/•M m [max] 

/ M m f(M m )dM m = p c n m , (5) 

•/M H i[min] 

where p c is the critical density at redshift z and SIhi is the contribution to the density parameter 
from the HI in DLAs at that redshift. Absorption studies (Lanzetta et al. 1995) have shown that 
in a flat cosmological model, in the redshift range z ^ 3.5, the total mass density contributed by 
the HI in DLAs can be fitted by the function 

n m (z) = 0.19 x Kr 3 /^ 1 exp(0.83z) . (6) 

We use this to fix the normalization constant A for the HI mass distribution function. 



2.2. Exponential disk (ED) 

Recent studies of the velocity structure of the DLAs show that the HI density across the cloud 
is consistent with a thick, rotating disk of exponential profile (Prochaska Sz Wolfe 1998) 

n(R,Z) =n exp - ^\ , (7) 

where n(R, Z) is the number density of HI atoms in the disk, R is the radial coordinate in the plane 
of the disk and Z is the coordinate along the thickness of the disk. Observations suggest that disks 
have Rd/Zd — 3, and are rotating with v c ~ 200-250 kms -1 . Excluding the effect of rotation, the 
gas in the disks has negligible velocity dispersion. Absorption studies do not fix R^ and they are 
only mildly dependent on the central column density A^hi(0) = 2 Z^nQ. In our analysis we assume 
that all the disks are observed face-on. As a consequence the column density distribution in each HI 
cloud is circularly symmetric around the center and falls off exponentially in the radial direction. 
For v = 320 MHz, the flux from a HI cloud is 

i\ / d. \ 2 



M HI ^ X1 O-M f *W U*^. (9) 



and the HI mass is 

A/Tttt = r v i n 10 1U„ l 

10 22 cm~ 2 y \10 kpc. 

We choose the central column density, Ahi(0) to have a fixed value in the range 5 x 10 21 cm~ 2 to 
2 x 10 22 cm -2 and allow R& to vary from 2 kpc to 10 kpc. For a fixed value of AT H i(0) = 10 22 cm -2 
this corresponds to a mass range Mni[min] = 2 x 10 9 M Q and Mhi [max] = 5 x 10 10 M Q . The 
treatment of the mass distribution function for the ED model has been carried out in exactly the 
same way as for the UD model, and we consider values in the range 0.5 to 1.0 for the slope a. 
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2.3. The flux distribution of HI clouds 

For both the UD and the ED model the HI mass distribution function can be converted to a 
flux distribution function using either Eqs. (3) and (4) or Eqs. (8) and (9). This requires assuming 
a value for AV which we take to be 200 kms" 1 . The GMRT channel at 320 MHz (z = 3.4) has 
a bandwidth of 16 MHz. This corresponds to a redshift range of 3.33 to 3.55. We use the flux 
distribution function to calculate the number of clouds per unit flux per arcsec 2 within the redshift 
range covered by the GMRT bandwidth. 

For the UD model where the value of tdla is fixed this gives us 



1.2 x 10 >Jy arcsec x 



dF u dVt 



. (10) 

iok PC ; yi^yj y ! 

For this model the different fluxes correspond to different values of the column density. For the 
ED model where the central column density 7V HI (0) is fixed, the number of clouds per unit flux per 
arcsec 2 is 

di ft A 1 O 

4.25 x 10 fxJy arcsec x 



dF v <Kl 

(In Eqs. (10) and (11) the normalization has a weak dependence on the parameter a, which is not 
shown for brevity) In this model the variation of flux corresponds to different values of the radial 
scale length For v = 320 MHz, the total number of clouds expected above a given flux in a 
1' x 1' region of the sky and within the GMRT bandwidth is shown in Figure 1 for the UD and ED 
models. 



3. Gravitational Lensing by clusters 

Here we give a brief summary of the salient features of gravitational lensing relevant to this 
work. A more thorough exposition can be found in Schneider et al. (1992). 

An object which is gravitationally lensed appears at a position which is different from its 
unlensed position. This information is encoded in the lens equation which relates the unperturbed 
position of the source to its perturbed position. If the angular coordinates of the source are j3 and 
the angular coordinates of the image are 0, then the lens equation is given by 



(3 = 6- W(0) 



(12) 
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The dimensionless relativistic lens potential ip satisfies the two-dimensional Poisson equation V 2 ip{6) = 
2k(9), where the convergence k{6) = £(#)/£ cr and is the two-dimensional surface mass den- 
sity of the lens, S cr = (c 2 / 'AitG)(d s / 'didi s ) being the critical density. The quantities d s , di, and di s 
are the angular-diameter distances from the observer to source, the observer to lens, and the lens 
to source respectively. Gravitational lensing does not effect the specific intensity along a light ray. 
The flux received by an observer is proportional to the solid angle subtended by the image at the 
observer; and since the solid angle of the image after lensing is, in general, different from that of 
the source, an observer can receive more (or less) flux in the lensed case than in the unlensed case. 

The shape and size of the image are related to those of the source by the transformation 
matrix M^ 1 = dPi/d9j. The infinitesimal area d 2 6 in the image plane is mapped to an area 
d 2 f3 = n~ l {6)d 2 6 in the source plane, where \x = det[My] is the magnification. The regions 
(curves) in the source plane where det[d Pi/dOj] = are called caustics. Magnification is infinite for 
point sources which are placed on the caustics. For finite sources the amplification remains finite, 
though the caustics still remain the points in the source plane where the magnification is large. 

For this paper we model clusters of galaxies as single component Pseudo Isothermal Elliptical 
Mass Distribution (PIEMD) after Kassiola & Kovner (1993). These are characterized by a core 
radius R c , one-dimensional velocity dispersion a, and ellipticity parameter e. In addition to these 
parameters, the lens redshift z\ and the source redshift z s completely specify the dimensionless 
potential ip. Real clusters can have substantial sub-structure and are best modeled as multi- 
component PIEMDs, but for simplicity we consider only single component lenses. 

For our case the source redshift, z s = v e f v Q —\ with v e = 1420 MHz, and v Q = {233, 320, 620} MHz. 
For our study we consider the following range of parameters: 1200 kms" 1 < a < 1800 kms -1 ; 
40 kpc < R c < 60 kpc; < e < 0.3, and 0.1 < z\ < 0.3. This choice of parameters is motivated by 
the observed properties of known clusters and their abundances (e.g. Mazure et al, 1996). 



The flux from a lensed HI cloud can be expressed in terms of the source coordinates (3 and the 
magnification fi((3) as: 



3.1. Lensing of HI clouds 



F, 




(13) 



This can be expressed in terms of the flux of the unlensed source 



F v [Lensed] = F u [Source] x /i , 



(14) 
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where fi is the magnification averaged over the angular extent of the source 

£(0, R) = J d 2 (3' ntf - (3') W(&, R) . (15) 

The normalized window function W(j3,R) represents the radial profile of the HI cloud and for the 
UD and ED models we use 



W((3,r DLA ) = _^9(|/3|-^M (UD) 



^DLA V ds J 



^Lexpf-^i 
2nR\ v \ R d / ds 



W((3,R d ) = —^expi—^) (ED), (16) 



where is the Heaviside step function. In our calculations we have evaluated p,(/3,R) by first 
calculating fx(J3) on a finely-spaced grid in the source plane. This is then convolved with the window 
function using the Fast Fourier Transform algorithm (Press et al. 1992) to obtain p,(J3,R). We 
have tested the convergence of this procedure with respect to variation of the grid size. In situations 
where there are multiple images we have added up the magnification of all the magnified images; 
because the typical image separation is less than the beam width of the GMRT and the observations 
which we are discussing here will not be able to resolve the multiple images. 

We use the procedure discussed above to calculate p,(/3, R) on a grid. This is used to calculate 
A(> jlo) (p,Q is any fiducial average magnification) which is the area (in arc-second 2 ) in the source 
plane for which the condition p,(j3, R) > £o is satisfied. For observations with a detection threshold 
of -Fjy[min], a HI cloud with unlensed flux F v will have a lensed flux greater than F^[min] for any 
position within the area A{> F u [mm]/F u ). Summing up the contributions from HI clouds with all 
possible values of the unlensed flux gives us Af(>F u [mm]) the total number of HI clouds that can 
be detected in the field of a cluster 

tf{>F„[mm]) = J dF u A (> F u [min]/F„) . (17) 

It should be pointed out that throughout the analysis we have assumed all the HI clouds to 
be oriented face-on. In reality the HI clouds will occur at all possible angles. In this case the 
mean magnification fi((3, R) will differ with the orientation of the cloud. However, the projected 
area is maximum when the cloud is face-on and the mean magnification will be lower relative to 
the situation where the cloud is located at the same position with an edge-on orientation. This 
means that our analysis gives an underestimate of the magnification and we may expect larger 
magnifications when all possible orientations are taken into account. 



4. Results 



We first discuss our results for observations centered at 320 MHz. The results for 620 and 
233 MHz are summarized at the end of this section. 
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We have calculated A(>p,) and jV(>i fl I/ [min]) for several cluster parameters with the disk 
parameters varying over the allowed range discussed earlier for the UD and ED models. The effect 
of varying the cluster parameters is to change A[>p) which in turn affects the number of detected 
HI clouds. The regions of high magnification (ji > 10) are restricted along the caustics and all 
the contribution to A(>jl) is from a narrow strip along the caustics. Although the length of the 
caustic curve increases as the ellipticity e is increased, this is accompanied by a fall in the values 
of the magnification. The value e = 0.1 is a good compromise between these two effects and the 
results are presented for this value. The e dependence of the results is rather weak in the range 
0.1 < e < 0.3. The values of A{>Ji) and jV(>i^[min]) decrease by around 10% as the value of e 
changes from 0.1 to 0.3. Amongst the cluster parameters, changing the velocity dispersion a has 
the strongest effect on A(>p) and this is shown in Figures 2 and 3. Varying the cluster redshift 
from 0.1 to 0.5 causes A(>p) to decrease by nearly 30%. 

Varying the radius of the HI clouds has a few distinct effects. First, an increase in the size of 
the HI clouds decreases the average magnification p,((3,R). The area A(>p) also decreases as the 
radius of the HI cloud is increased and this is clearly seen in Figures 2 and 3. The second effect is 
that a larger cloud radius gives a larger value of the unlensed flux (Eqs. ( 3) and (8)). The third 
effect is that as the radius of the HI clouds is increased, the number of clouds in the field decreases. 
This is seen in Figure 1 for the UD model. All these effects combine to determine how the number 
of detected HI clouds depends on the radius of the clouds. As is seen in Figure 4 for the UD model, 
an increase in the cloud radius generally increases the number of detected HI clouds. Only for 
large values of tdla do we find a decrease in the number of HI clouds causing a decrement in the 
number of detections. The results for the ED model shown in Figure (5) are qualitatively similar 
to those for the UD model. In this case the parameter that controls the number of clouds in a 
field is AT H i(0), and the number of clouds in a field falls as Ahi(0) is increased (Figure 1). As in 
the UD case, this generally leads to an increase in the number of detected clouds because a fall in 
the number of clouds is often over-compensated by an increase in the unlensed flux. An important 
parameter in both the cloud models is the spectral index, a (Eqs. (10) and (11))- The number of 
detected clouds decreases by 25% as a is increased from 0.5 to 1. 

For a detection threshold of 50 n Jy, for observations centered on a cluster with a = 1200 km s _1 , 
the expected number of detected clouds is in the range 0.1 to 0.5 for tdla > 10 kpc in the UD model 
and for ./Vhi(O) > 5 x 10 21 cm -2 in the ED model. For both the UD and the ED models the lower 
limit for detection approximately corresponds to the mass range 5 x 10 s M < Mm < 2.5 x 10 10 M Q . 

The number of detections will go up if the observations are centered on a cluster with a higher 
velocity dispersion, since M(> 50 //Jy) goes up by a factor of 2 to 3 if a = 1800 kms -1 instead of 
a = 1200 kms -1 (Figures 4 and 5). In addition it will be possible to detect the HI clouds even if 
they have lower HI masses. 

Future radio telescopes will reach sensitivities of around 2 /xJy. Eqs. (3) and (8) show that it 
may then be possible to detect the HI clouds without the aid of gravitational lensing. In the UD 
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model if tdla > 20 kpc then a substantial fraction of the HI clouds will be detected without any 
gravitational lensing. While a small fraction of the HI clouds may be detected without gravitational 
lensing if 10 kpc < tdla < 20 kpc, it will not be possible to observe the clouds without the aid of 
gravitational lensing if tdla < 10 kpc. Figure 6 shows that observations with a detection threshold 
of 2^Jy centered on a cluster with a = 1200 kms" 1 are expected to have a few detections even 
if the cloud radius is as small as tdla = 2 kpc. This corresponds to a mass range 2 x 10 7 M < 
M ffl < 10 9 M . 

The discussion until now has been restricted to 320 MHz. The main differences which arise at 
other frequencies are: (a) the flux from HI clouds changes (Eq. (2)), (b) the total number of clouds 
is different owing to changes in Ohi (Eq. (6)) and Az and, (c) the magnification of sources changes 
owing to change in the source redshift (§3). 

At z — 1.3 (610 MHz) observations indicate a ~ 2.5 and iVmfmax] ~ 10 21 cm -2 (Lanzetta et al. 
1995). Also, Ohi falls considerably from its value at z ~ 3.4 However, owing to obscuration of quasar 
light by the dust present in HI clouds (Fall & Pei 1995) these values may not be representative of 
the population of HI clouds and it is possible that the observed Ori is underestimated by a factor 
of 3 (Fall & Pei 1995) at this redshift . Given the uncertainty in the HI content we have used 
several values of spectral indices and Ori- Figure (7) shows the results for just one set of cluster 
parameters. 

The value of Qui at z ~ 5.10 (233 MHz) is not known. For the purposes of this paper we use 
the value given by Storrie-Lombardi et al. (1996) which is valid for z < 4.7. For v = 233 MHz the 
number of detectable clouds in a cluster field is shown in Figure (7). 

5. Discussion and Conclusions 

In this paper we have investigated the possibility of detecting the redshifted 21 cm emission 
from HI clouds at high redshift. Such observations are not possible with existing radio telescopes 
unless we observe the cloud through a cluster gravitational lens which magnifies the HI flux. How- 
ever, this method has the disadvantage that it will work only for clouds which lie very close to the 
caustics of cluster lenses. 

Given the lack of a clear picture about the nature of these objects we have used two simple 
models, namely the UD and the ED models (§2). The results are similar for both the models 
(Figures 4 and 5). Both the models have free parameters which effectively allow us the freedom 
of distributing the total HI available in the HI clouds at high redshift into either a large number 
of clouds with a small amount of HI each or a small number of clouds with large amounts of HI 
each. In the former scenario the probability of a cloud being located very close to the caustic of 
a gravitational lens and experiencing a large magnification is higher but this is offset by the fact 
that the unlensed flux of the HI clouds will be very small and a magnification of ~ 50 may not be 
sufficient to render the lensed flux above the threshold value for detection. The parameter range 
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where the probability of a detection is maximum changes with the threshold flux and this may lead 
to the possibility of using such observations to constrain the distribution of the HI masses of the 
clouds. This possibility has not been studied here in any detail. 

We find that for observations at 320 MHz with a detection threshold of 50 /uJy centered on 
a cluster with a = 1200 kms -1 , the chance of detecting a HI cloud is greater than 10% provided 
the minimum HI content of these clouds is in the mass range 5 x 10 s M Q < Mm < 2.5 x 10 10 M Q . 
The chances of detecting an HI cloud increases for a cluster with a higher velocity dispersion, and 
a single deep image in the field of a cluster such as Abell 1689 (Tyson et aZ.1990; Pierre et a/.1991; 
Smail et al.1991), which has a = 1989 kms -1 and zi = 0.196, might either reveal a cloud or put 
meaningful bounds on the mass range and mass spectral index a of the HI clouds. 

At 320 MHz the GMRT will reach a sensitivity of 50 /iJy with 100 hrs of integration at a 
frequency resolution of Au = 1.25 x 10 5 Hz corresponding to a velocity width Ay ~ 115 kms -1 . 
An alternative strategy would be to observe several clusters to a threshold flux of 100 /iJy. The 
latter strategy may be superior for a part of the parameter range shown in Figures 4, 5 and 7, i.e. 
the number of expected detections will be higher for the same total observation time. 

Future radio telescope Square Kilometer Array (SKA) will reach sensitivity ~ 2 //Jy at v ~ 
320 MHz for Ay ~ 200 kms -1 in an integration time of ~ 8 hrs 1 . In this case observations 
centered on clusters with a = 1200 kms -1 are expected to detect a few HI clouds even if the clouds 
have low HI masses in the range 2 x 10 7 M < M m < 10 9 M . 

Our results indicate that the probability of detecting a gravitationally lensed HI cloud at 233 
and 610 MHz is generally lower than the probability at 320 MHz. However, it should be borne in 
mind that there are larger uncertainties in our predictions at 233 and 610 MHz as compared to the 
predictions for 320 MHz. 

We conclude by noting that the strategy of carrying out deep radio observations centered on 
known cluster lenses will make it possible to detect the redshifted 21 cm emission from lensed HI 
clouds with either the GMRT or the SKA, depending on the HI content of these clouds. Such 
observations will shed new light on these objects whose HI content has been observed only in 
absorption to date. 

The authors would like to thank Jasjeet Bagla, Jayaram Chengalur, Divya Oberoi, and Somak 
Raychaudhury for useful discussions and an anonymous referee for helpful comments. T.D.S thanks 
CSIR for providing financial support. 



1 for details see http://www.nfra.nl/skai/science 
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Appendix A 

Here we present a derivation of Eq. (1) which gives the specific intensity of the redshifted HI 
emission from a cloud which has HI column density 7V HI , velocity dispersion AV and is at a redshift 
z. 

The HI contained in a part of the cloud which subtends a solid angle Af2 at the observer is 
Nmr\(z)AQ, where r\(z) is the angular diameter distance to the HI cloud. In the rest frame of 
the cloud the total luminosity of the HI emission is 

AL = A 21 h P u e fN m r 2 A (z) An , (18) 

and it will be distributed over the frequency interval Av e = AVv e /c. Here v e = 1420 MHz is the 
rest frame frequency of the HI emission, hp is the Planck constant, / is the fraction of HI atoms 
in the excited state and A 2 \ = 2.85 x 10~ 15 s~ l is the Einstein coefficient for the spontaneous 
emission. 

At the observer the radiation will be redshifted to a frequency v = v e j(\ + z) and it will be 
distributed over the frequency interval Au = Av e /(1 + z) The flux AF is given by 

AF = , (19) 

where rp(z) is the luminosity distance to the HI cloud. We use this to calculate the specific intensity 

AF 

= A^n ■ < 20 > 

Using Eqs. (18) and (19) and the relation (see e.g. Peebles 1993) 

ta(z) 
tl{z) 

we obtain 



(i+*r, (2i) 



T _ N m fA 2 i h P c 
v ~ 47T (1 + zf AV 1 > 

In a situation where the spin temperature T s is high (T s » hpu e /kp) we have / = 3/4, where 
hp denotes the Boltzmann constant. Equation (22) can be evaluated to give 

„ T 2/ n -\( N m \ /200 kms- 1 \ , , 

/, = 40 /xJy arcsec- 2 (l + z)' 3 ( ^ ™ _ 2 j f Ay j . (23) 
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Fig. 1. — This shows the number of HI clouds with flux greater than F v in a 1' x 1' patch of the 
sky for a bandwidth of 16 MHz centred at 320 MHz, and a mass spectral index a = 0.5. Left Panel 
is the UD model for three values of tdla (in kpc). Right Panel is the ED model for three values of 
iVffl(O). 
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Fig. 2. — The area in which the averaged amplification p, exceeds a given value \i is shown for 
several cluster parameters and clouds radii in the the UD model. In each panel, the clouds radius 
Hdla = {10,20,30,40} kpc corresponds to curves of decreasing maximum amplification. 
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Fig. 3. — The area in which the averaged amplification Jx exceeds a given value [i is shown for 
several cluster parameters and radius parameter Rd in the the ED model. In each panel, the radius 
parameter = {2, 8, 14} kpc, corresponding to curves of decreasing maximum amplification. 
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Fig. 4. — This shows the number of lensed HI clouds expected to be detected in the field of 
a cluster as a function of the cloud radius in the UD model with <x = 0.5, for the 4 cluster 
models shown in Figure 2. The three curves in each panel correspond to detection thresholds 
F^fmin] = {50,100,150} /i/Jy, with the number of expected detections decreasing for higher flux 
thresholds. 
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Fig. 5. — This shows the number of lensed HI clouds expected to be detected in the field of a 
cluster as a function of the central column density A^hi(O) in the ED model with a = 0.5. for 
the 4 cluster models shown in Figure 3. The three curves in each panel correspond to detection 
thresholds F^fmin] = {50,100,150} fiJy with the number of expected detections decreasing for 
higher flux thresholds. 
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Fig. 6. — This shows the number of lensed HI clouds expected to be detected in the field of a 
cluster as a function of the cloud radius, tdla, in the UD model with a = 0.5 for observations 
with Fjyfmin] = 2 /xJy and a bandwidth of 16 MHz centred at 320 MHz. The cluster parameters 
are shown in the figure. 
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Fig. 7. — This shows the number of lensed HI clouds expected to be detected in the field of a 
cluster with a = 1800 kms -1 , e = 0.1, z\ = 0.2 and R c = 45 kpc in the UD model. Left Panel is at 
610 MHz. The solid curve is for N m [max] = 10 21 cnT 2 , a = 2.3 and Q m given by Eq. (6). The 
dashed (dot-dashed) curve uses a = 2.3 (a = 0.5) with fini twice the prediction of eq. (6) and 
Nm[max) = 10 22 cm~ 2 Right Panel is at 233 MHz with a = 0.5 (solid line) and a = 1.5 (dashed 
line) 



